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ABSTRACT 


A  one-dimensional  mixed  layer  model  (Garwood,  1976,  1977)  is  used  to 
simulate  the  mixed  layer  depth  and  temperature  observed  during  POLEX,  a 
component  of  the  North  Pacific  Experiment.  POLEX  occurred  during  Janu- 
ary-February  1974  with  mid-ocean  observations  (35°N,  155°W)  of  tempera¬ 
ture  and  salinity  made  from  R/P  FLIP,  while  FLIP  was  under  free-dri  ft 
conditions.  The  results  of  this  simulation  show  that  apparent  horizon¬ 
tal  advection,  due  to  the  drift  of  R/P  FLIP,  was  important  in  the  heat 
and  salinity  budgets,  but  that  vertical  mixing  and  mixed  layer  depth 
changes  were  controlled  primarily  by  the  one-dimensional  response  of  the 
turbulent  kinetic  energy  budget  to  the  local  atmospheric  forcing.  Oc¬ 
casionally,  surface  salinity  flux,  due  to  large  precipitation  rates,  can 
significantly  alter  (on  the  order  of  tens  of  meters)  the  thickness  of 
the  ocean  surface  turbulent  boundary  layer,  with  a  demonstrated  decrease 
in  thickness  during  a  single  period  of  strong  precipitation. 
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I.  INTRODUCTION 


A.  PURPOSE  OF  STUDY 

The  objective  of  this  research  was  to  investigate  the  conclusions 
which  were  reached  by  Simpson  and  Paulson  [1979]  from  ocean  data  ob¬ 
tained  as  part  of  the  Pole  Experiment  (POLEX),  a  component  of  the  North 
Pacific  Experiment  (NORPAX).  POLEX  was  part  of  the  first  process  ori¬ 
ented  NORPAX  experiment.  The  experiment  was  named  POLE  to  indicate  that 
the  horizontal  extent  of  the  sampling  was  small  compared  to  the  largest 
scales  investigated  in  NORPAX.  Simpson  and  Paulson  concluded:  "One¬ 
dimensional  mixed-layer  deepening  models  failed  to  predict  the  mixed- 
layer  depth  and  temperatures  observed  during  POLE.  Horizontal  advec- 
tion,  as  evidenced  from  the  salinity  maximum  frequently  occurring  at  the 
bottom  of  the  mixed  layer  and  other  near-surface  changes  in  salinity  and 
temperature  not  associated  with  local  surface  forcing,  are  responsible 
for  the  failure".^ 

This  study  will  investigate  the  one-dimensionality  of  the  POLE 
temperature  and  salinity  data  and  determine  if  the  mixed  layer  depth 
changes  in  response  to  the  local  atmospheric  conditions  (one-dimensional) 
or  if  these  changes  are  caused  by  two  and/or  three  dimensional  advective 
processes.  To  accomplish  this  task,  a  non-stationary ,  one-dimensional 
bulk  model  of  the  mixed  layer  originally  proposed  by  Garwood  [1976,  1977] 


Simpson,  J.  J.  and  Paulson,  C.A. ,  1979:  Observation  of  Upper  Ocean 
Temperature  and  Salinity  Structure  During  POLE  Experiment,  J.  Phys. 
Qceanoqr. ,  9,  page  869. 
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is  expanded  by  including  a  salinity  budget.  A  comparison  of  the  re¬ 
sults  of  a  strictly  one-dimensional  model  of  the  upper  ocean  with  obser¬ 
vations  can  be  used  to  make  an  assessment  of  the  one-dimensionality  of 
the  different  variables  which  are  involved  in  the  upper  ocean  boundary 
layer  response  to  the  local  atmospheric  forcing. 

One  hypothesis  to  be  examined  is  that  horizontal  advection  may  be 
important  for  the  heat  and  salt  budgets,  and  at  the  same  tiVe  vertical 
mixing  and  mixed  layer  depth  may  be  controlled  primarily  by  one¬ 
dimensional  or  local  forcing  alone.  One  objective  of  this  study  will 
be  to  determine  how  responsive  are  the  variables  of  heat  and  salt  in 
the  upper  ocean  boundary  layer  to  one-dimensional  effects.  A  second 
objective  of  this  study  is  to  use  the  data  obtained  as  part  of  POLEX  to 
see  if  the  one-dimensional  mixed  layer  model  proposed  by  Garwood  will 
predict  the  observed  mixed  layer  depth,  even  if  advection  dominates  the 
local  changes  in  heat  and  salt  content  in  the  upper  ocean. 

The  data  used  in  this  report  contain  observations  made  from  R/P  FLIP 
while  under  free-dri  ft  conditions.  As  a  result  of  the  free-dri  ft  condi¬ 
tion  of  FLIP  during  the  data  collection  period,  advection  as  seen  by  R/P 
FLIP  will  be  a  combination  of  the  actual  advection  in  an  Eulerian  coor¬ 
dinate  system  for  the  area  and  advection  due  to  the  drift  of  FLIP. 
Henceforth,  the  term  advection  will  be  used  to  mean  a  combination  of  the 
actual  advection  term  and  apparent  horizontal  advection  due  to  the  drift 
of  FLIP. 

The  data  were  collected  during  the  period  28  January  74  through 
14  February  74.  During  that  time,  FLIP  occupied  a  station  approxi¬ 
mately  800  miles  north  of  the  Hawaiian  Island  Chain  under  free  drift 
conditions.  The  position  of  FLIP  ranged  from  35°39'N,  155°05'W  to 
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34°36'N,  1 55 °25 ' W.  The  experimental  site  is  hydrodynamical ly  complex, 

as  shown  in  Figure  1.  The  Subtropical  Front  is  known  to  meander  between 

32°  and  35°N  (Roden,  1974;  Barnett, 1976).  The  region  of  the  trade  winds 

northeast  of  Hawaii  has  air-sea  fluxes  of  latent  heat  in  excess  of  35 
-2  -1 

cal  cm  day  (Wyrtki,  1965).  The  Subtropical  water  mass  formed  in  this 
region  contrasts  markedly  with  the  less  saline  Eastern  North  Pacific 
Central  Water  characteristically  encountered  north  of  35°N.  The  Horse 
Latitudes  are  located  only  3°  of  latitude  to  the  south  of  the  observa¬ 
tion  area  and  the  North  Pacific  Current  is  expected  to  affect  the  gen¬ 
eral  hydrography  of  the  region. 

B.  MILITARY  IMPORTANCE 

Mixed  layer  depth  is  an  important  factor  when  dealing  with  sound 
propagation  in  the  ocean.  It  is  used  to  determine  such  sonar  detection 
variables  as  convergence  zone  distances,  source  level  depth,  and  cutoff 
frequency  for  surface  duct  propagation.  A  model  which  could  predict  the 
change  in  mixed  layer  depth  in  response  to  the  local  atmospheric  forcing 
could  be  incorporated ‘in  the  existing  acoustic  models,  for  example  the 
Fast  Asymptotic  Coherent  Transmission  (FACT)  and  the  Naval  Interim  Sur¬ 
face  Ship  Model  (NISSIM),  presently  being  used  at  Fleet  Numerical  Ocean¬ 
ography  Center  (FNOC).  A  one  dimensional  mixed  layer  model  would  help 
the  naval  tactician  to  directly  correlate  the  variability  of  the  ocean 
with  his  weapon  system  and  determine  the  most  efficient  deployment  of 
the  naval  Anti-Submarine  Warfare  (ASW)  assets. 

C.  LITERATURE  REVIEW 

The  vertical  fluxes  of  heat,  salt,  and  momentum  across  the  sea-air 
interface  are  the  sources  of  almost  all  oceanic  motions.  In  the  fully 
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turbulent  oceanic  mixed  layer  bounded  by  the  sea-air  interface  above  and 
by  the  dynamically  stable  water  mass  below,  the  vertical  fluxes  are 
large.  Below  the  mixed  layer,  they  are  usually  negligibly  small  so  that 
one  can  decouple  the  mixed  layer  from  the  underlying  stable,  quiescent 
water  mass.  This  homogeneity  is  the  root  of  the  term  "slab",  which  is 
often  used  to  describe  the  layer  as  depicted  in  Figures  2  and  3. 

When  dealing  with  sound  propagation  in  the  ocean,  the  term  mixed 
layer  depth  is  sometimes  taken  to  be  that  area  with  a  uniform  tempera¬ 
ture  structure,  but  this  is  not  always  the  case.  In  Figure  4  the 
bottom  of  the  mixed  ,_yer  depth  (Region  I)  appears  to  be  at  about  80 
meters,  and  below  there  is  a  strong  temperature  gradient.  However,  ex¬ 
amination  of  the  ot  profile  shows  that  the  mixed  layer  depth  is  about  30 
meters.  This  shallow  halocline  was  probably  caused  by  an  earlier  rain¬ 
fall  which  decreased  the  near-surface  salinity.  The  assumption  of  ver¬ 
tical  homogeneity  in  a  bulk  model,  because  of  the  problem  associated 
with  solving  for  the  interior  fluxes  of  buoyancy  and  momentum,  reduces 
to  the  need  to  know  only  the  surface  and  entrainment  fluxes.  However, 
only  small  vertical  gradients  in  these  mean  variables  may  be  associated 
with  large  turbulent  fluxes.  Therefore,  the  slab  assumption  should  not 
be  as  automatically  applied  to  the  turbulent  kinetic  energy  budget. 

The  immediate  local  reaction  of  this  mixed  layer  to  those  fluxes  re¬ 
sults  in  a  homogeneous  water  column,  i.e.,  vertical  uniformity  in  the 
mean  velocity  and  density.  This  adjustment  to  the  density  and  velocity 
structure  of  the  surface  layers  of  the  ocean  to  variable  fluxes  has  been 
the  subject  of  a. large  number  of  studies  since  Ekman's  [1905]  treatise 
where  he  originated  the  concept  of  a  depth  of  frictional  resistance  for 
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the  upper  section  of  a  wind  stressed  ocean.  He  further  suggested  that 
the  wind-driven  current  will  have  a  similar  depth,  the  Ekman  depth. 

Much  of  the  one-dimensional  theory  for  the  ocean  surface  turbulent 
boundary  layer  or  mixed  layer  seems  to  be  dependent  upon  the  validity 
of  two  basic  hypotheses.  The  first  of  these  is  that  vertical  mixing 
within  the  turbulent  boundary  layer  and  entrainment  mixing  at  its  base 
occur  in  response  to  local  atmospheric  forcing,  i.e.,  the  surface  wind 
stress  and  net  heat  radiated  or  fluxed  across  the  sea-air  interface.  The 
second  hypothesis  is  that  the  mechanical  energy  budget  is  the  key  to  the 
understanding  and  prediction  of  the  mixed  layer  dynamics. 

Most  observations  show  that  local  atmospheric  forcing  such  as  an 
increase  in  wind  speed  increases  the  mixed  layer  depth  and  low  wind  speed 
and  high  solar  radiation  decreases  mixed  layer  depth.  This  is  less  true 
with  increased  depth  below  the  mixed  layer.  Observations  by  Pollard  and 
Tarbel  [1975]  showed  that  horizontal  coherence  between  the  fluctuating 
wind-driven  currents  near  the  surface  was  large  and  consistent  with  the 
local  amplitude  and  phase  differences  of  the  wind  history.  With  in¬ 
creased  depth,  the  horizontal  coherence  decreased  monotonically  as  did 
the  vertical  coherence  of  the  deeper  current  measurements  with  those 
near  the  surface.  Kroll  [1975]  pointed  out  that  the  time  scale  of  these 
strongly  inertial  currents  is  on  the  order  of  a  couple  of  days,  depend¬ 
ing  upon  the  stability  of  the  underlying  water  column  as  well  as  the 
difference  between  the  driving  frequency  and  the  local  inertial  fre¬ 
quency.  The  turbulent  motions  within  the  mixed  layer  are  expected  to  be 
even  more  responsive  to  the  local  atmospheric  conditions  because  of  the 
small  dissipation  time  scale  for  the  turbulent  kinetic  energy — minutes 
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to  hours,  which  depends  upon  the  layer  depth  and  turbulence  intensity. 
Both  the  mean  and  turbulent  kinetic  energies  in  the  mixed  layer  are  ex¬ 
pected  to  reflect  the  local  boundary  conditions  for  wind  stress  and  sur¬ 
face  buoyancy  flux.  Measurements  of  turbulent  velocity  fluctuations, 
which  were  made  by  Powell  et  al  [1975]  and  Jones  and  Kenney  [1977],  show 
that  the  root-mean-square  turbulent  velocity  does  scale  with  the  water 
surface  friction  velocity,  u*.  Camp  and  Elsberry  [l 978]  and  Elsberry 
and  Camp  [1978]  demonstrated  the  role  of  strong  synoptic  scale  wind 
events  for  deepening.  Garwood  and  Hal  pern  [j  976]  and  Elsberry  and  Raney 
[1978]  indicated  the  importance  of  diurnal  heating,  especially  at  times 
between  storms.  Hindcasting  simulations  by  Thompson  [1976],  Mellor  and 
Durbin  [1975]  and  Camp  [l 9 7 6]  support  the  importance  of  surface  boundary 
fluxes. 

Kraus  and  Turner  |l967j  improved  the  one-dimensional  model  originally 
developed  by  Kraus  and  R ooth  [196]]  .  This  was  accomplished  by  consider¬ 
ing  the  turbulent  kinetic  energy  budget  utilizing  the  heat  equation  and 
a  mechanical  energy  equation.  Because  the  frictional  generation  of  heat 
is  negligible,  the  vertically  integrated  heat  equation  provides  a  rela¬ 
tionship  for  the  conservation  of  potential  energy.  They  parameterized 
the  mechanical  production  rate  in  terms  of  the  friction  velocity,  but 
neglected  the  viscous  rate  of  dissipation  and  the  effects  of  entrainment 
shear  production. 

The  application  of  mechanical  energy  principles  to  the  problem  of 
mixed  layer  dynamics  has  caused  some  disagreement.  An  important  contro¬ 
versy  in  mixed  layer  modeling  has  been  centered  around  whether  mean  kin¬ 
etic  energy  or  turbulent  kinetic  energy  leads  to  the  downward  entrainment 
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buoyancy  flux  at  the  base  of  the  turbulent  boundary  layer.  The  mean 
kinetic  energy  advocates,  including  Thompson  0 976)  ,  Pollard  0977]  and 
Price  0977] ,  argue  that  this  mean  flow  energy  is  directly  converted  to 
potential  energy  by  a  mean  flow  instability  that  causes  the  bulk 
Richardson  number  to  be  less  than  some  critical  constant  of  the  order 
one.  Turbulent  kinetic  energy  advocates  do  not  disagree  with  the  pos¬ 
sibility  of  this  mechanism,  but  they  argue  it  is  not  the  dominant  mech¬ 
anism  because  turbulence-generated  instabilities  will  usually  suffi¬ 
ciently  erode  the  interface  to  preclude  the  occurrence  of  the  mean  flow 
instability.  Pe  Szoeke  and  Rhines  0975)  and  Yun  0978]  demonstrated 
that  wind-shear  production  of  turbulent  energy  is  ultimately  the  most 
important  source  of  energy  for  deepening  the  mixed  layer,  if  the  surface 
buoyancy  flux  is  neglected.  Even  whin  the  mean  flow  energy  is  important, 
they  find  that  the  rate  of  conversion  of  this  energy  to  potential  energy 
via  entrainment  shear  production  is. controlled  by  the  constraints  of  the 
turbulent  kinetic  energy  budget.  The  Kraus-Turner  turbulent  energy  model 
has  not  proven  to  be  entirely  satisfactory,  and  a  number  of  changes  have 
been  suggested.  Stevenson  [_1 979]  examined  a  number  of  alternative  para- 
meterizations  for  dissipation.  Zeman  and  Tennekes  0975) ,  Elsberry  et  al 
M .  Alexander  and  Kim  0976]  and  Kim  0  976]  added  depth-dependent 
dissipation  terms  that  helped  to  reduce  unrealistic  deepening- 

Garwood  0976,  1977)  developed  an  ocean  mixed-layer  model  using  the 
Navier-Stokes  equation  of  motion  with  the  geostrophic  component  eliminated, 
the  continuity  equation  in  incompressible  water,  the  heat  equation  from 
the  first  law  of  thermodynamics,  the  conservation  of  salt  equation,  and 
a  linearized  equation  of  state.  The  entrainment  hypothesis  depends  upon 
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the  relative  distribution  of  turbulent  energy  between  horizontal  and 
vertical  components  and  is  offered  as  a  plausible  mechanism  for  govern¬ 
ing  both  entrainment  and  layer  retreat. 

Garwood  suggested  that  planetary  rotation  influences  the  dissipation 
of  turbulence  for  deeper  mixed  layers  and  enables  a  cyclical  steady  state 
to  exist  on  an  annual  basis.  Furthermore,  the  rate  of  entrainment  for  the 
stable  regime  can  not  be  a  simple  linear  extrapolation  of  the  unstable 
situations.  Unlike  the  atomospheric  boundary  case,  most  of  the  solar 
radiation  does  not  penetrate  the  layer.  Therefore,  downward  turbulent 
heat  flux  in  the  oceanic  boundary  layer  is  as  important  as  the  upward 
flux  during  the  course  of  both  diurnal  and  annual  cycles.  The  non¬ 
linearity  of  the  interface  entrainment  tendency  parameter,  which  is 
greatest  for  stable  surface  boundary  conditions,  results  in  a  modulation 
of  the  long-term  trend  of  the  mixed-layer  depth  by  the  diurnal  component 
of  surface  heat  flux.  In  this  model,  buoyant  production  is  somewhat 
more  efficient  than  shear  production  as  a  source  of  energy  for  vertical 
mixing  because  of  its  unique  effect  on  the  vertical  component  of  the 
turbulent  velocity.  The  buoyancy  equation  is  generated  from  the  heat 
and  salt  equations  together  with  an  equation  of  state, 

P-  =  pQ  [l  -  a  (8  -  0q)  +  S  (s  -  SQ)j  7  (1) 

and  the  definition  for  buoyancy, 

b  =  9  (p0  -  p)/p0  (2) 

In  equations  1  and  2  9  is  temperature,  S  salinity  and  p  density 
while  a  and  8  are  the  expansion  coefficients  for  heat  and  salt,  respec¬ 
tively,  and  g  is  gravity.  The  tilde  represents  the  total  instantaneous 
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value  and  the  subscript  zero  denotes  a  representati ve  but  arbitrary  con- 

%  % 

stant  value.  The  generalization  of  using  b  rather  than  0  will  cast  the 
model  equations  in  a  form  equally  applicable  to  those  situations  where 
evaporation  and  precipitation  contribute  significantly  to  the  surface 
buoyancy  flux  and  the  structure  of  the  evolving  pycnocline.  The  buoyancy 
equation  also  has  a  more  obvious  and  direct  role  in  the  mechanical  energy 
budget  as  shown  in  Figure  5. 

Miller  [j 976^1  included  salinity  in  the  mixed-layer  model  of  Kraus  and 
Turner  ['967]  .  He  concluded  that  salinity  is  important  in  determining 
the  density  structure  when  applying  a  mixed-layer  model  to  the  entire 
ocean,  but  the  relative  importance  of  salinity  effect  on  the  short  term 
evolution  of  the  density  profile  may  not  be  significant.  Paulus  [l978] 
included  salinity  in  the  mixed  layer  model  of  Camp  {l976].  He  concluded 
that  the  effects  of  salinity  structure  were  noticeable  only  in  deepening 
regimes. 
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II.  DATA  ANALYSIS 

A.  METHODS 

The  upper  ocean  temperature  and  salinity  data  and  the  meterological 
data  during  the  NORPAX  POLE  experiment  were  obtained  from  Simpson  and 
Paulson  p 9-77^ .  The  parameters  required  for  surface  boundary  condition 
computations,  which  are  used  in  Garwood's  model,  include  wind  speed,  wind 
direction,  cloud  cover,  sea  surface  temperature,  air  temperature  (dry 
bulb),  dew  point  and  rain  code.  These  parameters  were  recorded  at  ir¬ 
regular  time  intervals.  A  linear  interpolation  technique  was  used  to 
obtain  values  for  these  input  parameters  at  regular  time  intervals.  The 
resulting  interpolated  data  were  compared  with  the  original  data.  Fig¬ 
ures  6  through  12  show  that  no  noticeable  error  was  introduced  into  the 
parameters  by  the  use  of  a  linear  interpolation. 

Temperature  and  salinity  profiles  as  in  Figures  3  and  4  were  digi¬ 
tized  by  using  a  Houston  Instrument  Hi  Pad  Digitizer  and  a  Tektronix 
4052  computer.  The  digitized  profiles  of  temperature  and  salinity  were 
linearly  interpolated  into  hourly  intervals.  This  provided  for  both 
initial  conditions  in  the  one-dimensional  model  and  for  fields  to  be 
used  to  compute  observed  mixed  layer  depth,  h0(t),  and  to  compare  with 
subsequent  model  computations  of  T(z,t)  and  S(z,t). 

B.  SURFACE  FLUXES 

The  mixed  layer  model  uses  the  total  surface  heat  flux  ( Qy ( 0 ) ) ,  in¬ 
cident  solar  radiation,  wind  speed  at  a  height  d,  rate  of  evaporation 
(E),  and  rate  of  precipitation  (P)  to  calculate  surface  fluxes  of 
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buoyancy  (heat  and  salt)  and  momentum.  The  friction  velocity,  in  air. 


is  calculated  using  the  following  formulas: 


Ts  =  pacdua 


u*  =  ^Ts/pa^ 


1/2 


(3) 

(4) 


where  u.  is  the  mean  wind  speet  (m/sec), 

a 


-3, 


cd  is  the  non-dimensional  drag  coefficient  (1.4x10  ), 

3  3 

Pa  is  the  density  of  air  (1.25x10  gm/cm  ),  and 

o 

ts  is  the  surface  stress  (dynes/cm  ). 

The  turbulent  fluxes  of  latent  heat  (Qg)  and  sensible  heat  (Qh)  are 
estimated  using  the  following  bulk  aerodynamic  formulas: 

(5) 


%  •  cd<-98Es  -  Ea>“a 
«h  ’  cd<Ts  *  Ta>“a 


(6) 


The  net  back  radiation  (Qfa)  is  estimated  from  the  following  empirical 
formula  reported  by  Husby  and  Seckel  [1975]: 


Qb  =  1.14  x  10'7(273.16  +  Ts)4  (.39 


■  05Ey2)  (1.  -  .6C2)  (7) 

a 


where  E$  is  the  saturated  vapor  pressure  of  the  marine  air  directly  in 
contact  with  the  sea  surface  (.98  corrects  for  salt  defects), 

Ta  is  the  air  temperature  (degrees  Centigrade), 

E,  is  the  vapor  pressure  of  air  at  approximately  10  m  based  on  dew 

a 

point  temperature, 

T$  is  the  sea-surface  temperature  (degrees  Centigrade),  and 
C  is  the  fractional  cloud  cover. 

The  upward  heat  flux  is  then: 

Qa  -  Qa  +  Q,  +  Qk  (8) 


..  v 
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The  solar  Insolation,  Qs>  was  estimated  by: 

Qs  »  (1.  -  a*b)  (1.  -  . 66C3 ) Qq  (9) 

The  constants  a  and  b  are  adapted  from  Tabata  |_1 964]  and  the  cubic 
cloud  cover  correction  from  Laevast  Jj"96o| .  The  coefficient  C  is  the 
fractional  cloud  cover  and  <*  is  the  mid-day  elevation  angle  of  the  sun. 

The  clear  sky  radiation,  Qq,  is  given  by 

Qq=Aq+A^  cos4>+B-j  si  n<}>+A2COs24>+B2S  i  n2<{>  (10) 

from  the  formula  developed  by  Seckel  and  Beauday  |j973)  .  The  coeffi¬ 
cients  (Aq,  etc.)  were  calculated  by  harmonic  representation  of  the  values 
presented  in  the  Smithsonian  Meteorological  Tables  (List,  1958)  with 
4>=(2tt/365) (t-21 )  where  t  is  the  Julian  day  of  the  year.  Turbulent 
temperature  flux, 

(ID 


TV 


(0) 


-S2L 

PCP 


turbulent  velocity  flux, 
u'w' 


(0) 


•  U  * 


turbulent  salinity  flux. 


S'W(0)  =  <P-E)S(0) 


and  turbulent  buoyancy  flux, 

0TV(Q)  "  SS  ’ W *  ( 0 )] 


b'w' 


02) 


(13) 


(0)  =  9  |_ai  w  (0) 

where  a  is  the  coefficient  of  thermal  expansion, 

8  is  the  coefficient  of  contraction  of  sea  water  due  to  salinity 
variations,  and 

is  the  surface  salinity  value. 


(14) 
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Figures  13,  14,  15  and  16  are  graphs  which  show  the  hourly  values  of 
wind  stress,  t,  total  heat  flux  at  the  surface,  Qy,  salinity  flux,  S'w' , 
and  buoyancy  flux,  b'w',  respectively,  as  calculated  by  the  one-dimensional 
model.  These  figures  can  be  compared  with  Figures  17  and  18  which  show  the 
interpolated  hourly  values  of  wind  stress,  x,  total  heat  flux  at  the  sur¬ 
face,  Qy,  buoyancy  flux,  b'w',  and  various  other  components  which  were  re¬ 
ported  by  Simpson  and  Paulson  D9773  .  The  values  of  these  fluxes  cal¬ 
culated  here  are  nearly  identical  to  the  values  reported  by  Simpson  and 
Paulson.  In  Simpson  and  Paulson's  graph  on  buoyancy  flux  (Figure  17),  a 
positive  buoyancy  flux  results  when  Qngt  >  0.  In  the  one-dimensional 
model  (Figure  16)  a  positive  buoyancy  flux  results  when  Qnet  <  0  and  E  >  P. 

1 .  Surface  Wind  Stress 

The  wind  stress  calculated  using  the  bulk  approximation  with  a 

-3  -2 

drag  coefficient  of  1.4x10  has  a  mean  value  of  about  0.66  dynes  cm 

for  the  experimental  period.  Figure  13  shows  the  output  for  the  surface 

wind  stress  as  calculated  by  the  one-dimensional  model.  Analysis  of 

Figure  13  shows  that  for  periods  up  to  as  long  as  4  days  the  wind  stress 

_2 

is  considerably  less  than  1  dyne  cm  ,  and  it  is  evident  that  only  two 

periods  of  relatively  high  winds  were  encountered  during  the  experiment. 

One  period  of  high  wind  near  the  middle  of  the  experiment  (time  =  200-260 

-2 

hours)  with  a  maximum  wind  stress  of  2.5  dynes  cm  ,  and  the  other  at  the 

_? 

end  of  the  experiment  with  a  maximum  wind  stress  of  3.3  dynes  cm  . 

2.  Surface  Heat  Flux 

Comparison  of  the  total  heat  flux  with  net  all-wave  flux  (Figures 
17  and  19)  shows  that  the  net  all-wave  flux  is  generally  the  dominant 
term  in  the  surface  heat  balance.  However,  a  two-day  period  of  enhanced 


latent  heat  flux,  due  largely  to  a  significant  increase  in  wind  speed 
midway  through  the  experiment  (time  =  250  hours),  is  evident  when  com¬ 
paring  Figures  13  and  14. 

Simpson's  [1977]  calculation  of  the  total  heat  budget  for  the 

_2 

experiment  suggested  that  the  ocean  gained  201  cal  cm  for  the  time  per¬ 
iod  from  3  Feb  74  -  12  Feb  74.  Simpson's  heat  budgets  calculations  indi¬ 
cated  that  ocean  gained  heat  for  6  days,  was  in  near  thermal  equilibrium 
with  the  atmosphere  for  2  days  and  lost  heat  to  the  atmosphere  for  2-day 
period.  These  results,  coupled  with  the  low  wind  stress,  imply  that  the 
upper  dynamics  might  have  been  dominated  by  net  surface  heating  for  a 
significant  part  of  the  experiment. 

3.  Salinity  Flux 

Figure  15  shows  the  salinity  term  as  computed  by  the  one¬ 
dimensional  model  (equation  13)  without  including  the  precipitation  term 
in  the  computation.  The  precipitation  term  was  omitted  in  the  calculations 
because  no  reasonable  correlation  could  be  made  between  a  rain  code  given 
in  the  experiment  data  report  (Simpson  and  Paulson  1977)  and  an  amount  of 
precipitation.  The  results  of  model  runs  which  included  precipitation 
(estimated  by  the  rain  code)  were  compared  to  model  runs  which  neglected 
precipitation  altogether.  These  model  comparison  runs  showed  no  signifi¬ 
cant  difference  except  during  one  period  when  precipitation  was  much 
larger  than  evaporation.  Simpson  and  Paulson  0979^  observed  that  over¬ 
all  evaporation  and  precipitation  were  in  near  equilibrium,  with  precip- 
itation  exceeding  evaporation  by  23  mg  cm  . 

4.  Buoyancy  Flux 

Stratification  in  the  ocean  exerts  a  strong  dynamical  influence 
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on  mixing  since  a  stably  stratified  fluid  requires  work  against  buoyancy 
forces  if  mixing  is  to  occur.  For  a  stably  stratified  fluid,  the  density 
increases  with  depth  and  the  buoyancy  flux  is  directed  downward,  while 
the  converse  is  true  for  an  unstable,  stratified  fluid.  In  the  one¬ 
dimensional  model,  surface  buoyancy  flux,  b'w1 ^ ,  is  determined  by  the 
air-sea  transfers  (equation  14).  Figure  16  shows  the  model  computed 
buoyancy  flux  with  Qngt  <  0  and  E  >  P  resulting  in  a  positive  buoyancy 
flux.  The  net  heat  flux,  Figure  14,  dominates  the  surface  values  of 
buoyancy  flux  during  the  experiment.  During  daytime  periods,  negative 
buoyancy  flux  results  from  the  domination  of  the  solar  heating  at  the 
surface,  while  during  the  night  periods  a  positive  buoyancy  flux  is 
produced  by  the  combined  effect  of  net  long-wave  radiation  and  the  up¬ 
ward  turbulent  fluxes  of  heat  and  moisture. 


C.  ADVECTION 


The  determination  of  the  importance  of  advection  to  the  heat  and 
salt  budget  is  necessary  before  an  assessment  can  be  made  on  how  one- 
dimensional  is  the  response  of  the  variables  of  heat  and  salt  are  in 
the  upper  boundary  layer.  The  velocity  vector  V  of  the  advection  term 
(V‘(VT,VS)),  is  separated  into  two  components  in  order  to  determine 
which  component  dominates  the  advection  term. 


Heat  and  Salinity  Budgets 


Neglecting  internal  sources  and  sinks,  the  equation  for  conser 


vation  of  heat  or  salt  is: 


(15) 
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where  T  is  heat,  S  is  salinity,  t  is  time,  d  is  depth  and  V  is  the  velo¬ 
city  vector  in  a  coordinate  system  moving  with  the  drifting  platform. 

The  meanings  of  these  terms  are: 

(a)  Heat  or  salinity  change  with  time, 

(b)  Heat  or  salinity  change  due  to  advection  in  the  moving  coordinate 
system,  assuming  that  when  averaged  over  a  day  such  features  such 
as  internal  waves  would  not  contribute  to  advection, 

(c)  Vertical  heat  or  saltflux  divergence,  and 

(d)  and  (e)  Horizontal  heat  or  salinity  flux  divergence. 

Terms  (d)  and  (e)  are  assumed  to  be  small  because  they  can  not  be 
evaluated.  Therefore,  equation  (15)  can  be  rewritten  in  terms  of  the 
variable  temperature: 


<• 

fi+t-5T 

If 

d 

dz  =  -  (T'w'(0))  +  T'w'(d)  (16) 

(f)  (9) 

Term  (f)  is  equal  to  the  surface  heat  flux  Qy/pcp,  where  Qy  is  the  net 

O  _  I 

surface  heat  (Qj=Qs-Qt;)-Qe-Qh+Qc) ,  cal  cm  sec  ,  p  is  the  density  of  sea 

-3  -1  -1 

water,  1  gm  cm  ,  and  cp  is  the  specific  heat  of  sea  water,  1  cal  gm  C  . 

Term  (g)  is  assumed  to  be  negligible  if  d  is  well  below  the  mixed  layer, 

i.e.,  d  >  hm  .  The  vertical  integrated  mean  advection  of  temperature  is 

iTlaX 


therefore. 


(17) 
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The  same  procedure  used  to  derive  (17)  may  be  followed  to  find  the 


average  advective  change  in  salinity.  The  surface  salinity  flux  is 
given  in  equation  (13).  Thus  the  vertically  integrated  mean  advection 
of  salinity  is: 


The  magnitude  of  actual  advection  (in  an  Eulerian  coordinate  system) 

can  be  checked  by  estimating  V  and  VT  or  VS.  The  term  V‘(VT,  VS),  is  the 

advection  quantity  averaged  over  a  time  interval.  At.  The  velocity  vector 

V  can  be  separated  into  two  components, 

V  =  Vc  +  Vp  (19) 

where  V  is  the  velocity  vector  due  to  the  ocean  currents  and  V„  is  the 
c  p 

velocity  vector  due  to  the  drift  of  R/P  FLIP.  The  velocity  vector  due  to 
the  ocean  currents  is  itself  composed  of  two  components: 


(20) 


where  V  is 
9 


the  velocity  of  the  current  due  to  the  quasi-geostrophic  flow 


and  is  assumed  to  be  parallel  to  the  surface  temperature  contour.  There¬ 


fore  V^-VT  is  assumed  to  be  negligible.  The  order  of  magnitude  of  Vg, 
the  velocity  current  due  to  Ekman  flow,  is  a  function  of  the  wind  stress 


and  can  be  estimated  by, 

|Ve|  =  pg  =  1.57  cm  sec”^  =  0.0565  km  hour-^  (21) 

phf 

_  _2 

where  x  is  the  average  wind  stress,  0.66  dyn  cm  ,  p  is  the  density  of 
w 

sea  water,  h  is  average  mixed  layer  depth,  5000  cm  and  f  is  the  Coriolis 
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parameter,  2tDSin<}>,  where  w  is  the  angular  velocity  of  the  earth  and  <p 
is  the  latitude. 

The  velocity  calculated  due  to  Ekman  flow  is  insignificant  when  com¬ 
pared  with  the  drift  rate  velocity  calculated  for  R/P  FLIP  (Vp=0.3  km  hour"1) 
from  28  Jan  1974  to  14  Feb  1974.  The  calculated  Ekman  velocity  is  not  suf¬ 
ficient  enough  to  explain  the  drift  of  POLE.  Hence,  R/P  FLIP  was  by  no 
means  flowing  with  the  same  mixed  layer  water  mass  throughout  the  experi¬ 
ment.  Therefore,  the  velocity  vector  V*  will  depend  on  the  drift  rate  of 
R/P  FLIP,  V  =  V  .  What  Simpson  and  Paulson  0979*j  call  advection  was  not 
actual  advection  for  the  area  but  apparent  advection  due  to  the  drift  of 
R/P  FLIP. 

2.  Calculation  of  the  Rates  of  Advection 

The  distance  R/P  FLIP  traveled  from  0700  hours  28  Jan  1974  to 
0700  hours  14  Feb  1974  (408  hours),  was  calculated  to  be  121  km,  using 
a  rhumbline  distance  calculation  and  a  position  of  35°36'N,  155°05'W  to 
34°36N,  1 55°25 ' W.  Vertical  profiles  of  temperature  and  salinity  were 
taken  from  R/P  FLIP  from  0000  hours  30  Jan  1974  through  0700  hours  Feb 
1974  (total  of  367  hours)  which  covered  a  distance  of  approximately 
109  km.  Figures  19  and  20  show  the  observed  changes  in  heat  and  salt 
fluxes.  The  amounts  of  advection  in  a  coordinate  system  moving  hori¬ 
zontally  with  R/P  FLIP  were  then  estimated  by  subtracting  the  accumu¬ 
lated  surface  heat  and  salt  fluxes  from  the  changes  in  heat  and  salt 
content,  respectively.  Figures  21  and  22  show  the  results  of  this  cal¬ 
culation.  Least  square  fits  were  performed  on  the  results,  giving  a 
-2  -1 

slope  of  2.0  cal  cm  hour  for  heat  content  gain  and  a  slope  of 
0.45  °/00  cm  hour"1  for  salinity  gain.  These  values  can  be  used  to 
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compute  the  average  horizontal  temperature  and  salinity  gradients  by 
multiplying  by  the  total  hours  of  observation  (367  hours)  and  dividing 
by  the  average  mixed  layer  depth  (approximately  40m)  and  the  distance 
FLIP  traveled  during  the  observational  period  (109  km), 

|VhT|  *  0.0017  C  km'1  (22) 

jVHS|  *  0.0004  °/00  km'1  (23) 

where  |V^T|  is  the  horizontal  temperature  gradient  and  |V^S|  is  the  hori¬ 
zontal  salinity  gradient. 

For  comparison  purposes,  gradients  were  obtained  from  a  NORPAC 
Atlas  |j96oj  by  dividing  the  atlas  temperature  change  (AT  =  0.2  °C)  and 
the  atlas  salinity  change  (AS  =  0.05  ^/00)  between  FLIP  starting  and 
final  positions  (35°39'N,  155°05'W  to  34°36'N,  155°25'W)  by  the  total  dis¬ 
tance  (121  km)  which  FLIP  traveled, 

|VhT|a  =  0.0017°C  km'1 
|VHT|A  =  0.0004  °/00  km'1 

where  |V^T|A  and  | V^S ] A  are  the  horizontal  temperature  and  salinity 
gradients  that  were  calculated  from  the  NORPAC  Atlas.  Although  these 
values  agree  with  the  earlier  calculation,  equations  (22)  and  (23),  they 
were  only  obtained  to  show  agreement  in  order  of  magnitude  with  recorded 
atlas  values.  These  results  further  illustrate  the  fact  that  what 
Simpson  and  Paulson  Jj979j  call  advection  was  probably  the  result  of  the 
drift  of  R/P  FLIP. 

Figures  19  and  21  show  that  the  temperature  field,  which  is  a 
function  of  both  time  and  depth,  has  large  fluctuations  which  can  only 
be  explained  by  advection,  or  the  drift  of  FLIP  which  is  a  non  one¬ 
dimensional  effect.  These  large  fluctuations  dominate  the  heat  budget 
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on  the  time  scales  shorter  than  100  hours,  but  over  a  longer  time  scale, 
surface  flux  is  equally  important  (one-dimensional  effects).  Figure  23 
is  a  plot  of  observed  sea  surface  temperature  and  a  plot  of  the  model  - 
computed  mixed  layer  temperature,  with  the  model -computed  temperature 
slightly  offset  on  the  temperature  scale  for  comparison.  Comparing  the 
two  plots,  a  one-dimensional  component  can  be  detected  in  the  observed 
sea  surface  temperature  time  series  but  the  stronger  component  due  to 
the  drift  of  FLIP  is  evident.  Therefore,  this  apparent  advection  (as 
seen  by  FLIP)  of  temperature  is  important  on  all  time  scales,  but  domi¬ 
nates  the  short  time  scale.  Figures  20  and  22  show  that  advection  of 
salinity  dominates  the  salinity  budget  over  both  the  short  and  long  time 
scales  up  to  at  least  400  hours,  which  was  the  period  of  observations. 

D.  MIXED  LAYER  DEPTH 

1 .  Observed  Mixed  Layer  Depth 

The  observed  mixed  layer  depth,  shown  in  Figure  17  as  used  by 
Simpson  and  Paulson  jj979j,  is  the  shallowest  depth  at  which  the  den¬ 
sity  is  not  more  than  0.02  ot  units  greater  than  the  density  at  a  depth 
of  5  m.  However,  salinity  and  temperature  changes  can  occur  above 
z  =  5m.  In  such  cases,  the  above  definition  will  result  in  a  misleading 
computed  time  series  for  mixed  layer  depth.  A  more  practical  definition 
for  observed  mixed  layer  depth,  which  would  be  more  consistent  with  the 
precision  of  the  available  P0LEX  observational  data,  would  be  the  shal¬ 
lowest  depth  at  which  the  observed  value  of  at  is  0.02  at  units  greater 
than  the  observed  surface  value.  For  comparison  with  the  Garwood  model, 
the  observed  mixed  layer  depth  (hg)  will  be  defined  according  to  this 
latter  definition. 
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Contours  of  temperature  for  the  observation  period  are  shown  in 
Figures  24  and  25.  Figure  24,  which  is  based  on  a  contour  interval  of 
0.5  C,  shows  a  mixed  layer  depth  as  varying  from  near  the  surface  to  ap¬ 
proximately  90  meters.  Figure  25,  which  has  a  contour  interval  of  0.05  C 

2 

after  having  been  smoothed  twice,  shows  a  mixed  layer  from  near  the  sur¬ 
face  to  approximately  80  meters.  This  illustrates  the  importance  of 
resolution  when  determining  mixed  layer  depth  from  temperature  contours. 

A  higher  resolution  frequently  results  in  a  shallower  mixed  layer  depth. 
Therefore,  a  temperature  contour  interval  of  0.5  C  is  not  sufficient 
resolution  to  accurately  determine  mixed  layer  depth  from  the  contour 
plot.  Although  density  is  a  function  of  both  temperature  and  salinity, 
during  this  time  series  the  salinity  influence  was  not  significant,  so 
that  the  density  was  mainly  dependent  upon  temperature.  A  mixed  layer 
depth  definition  could  then  be  defined  as  the  shallowest  depth  at  which 
the  temperature  is  not  more  than  0.1  C  less  than  the  temperature  at  the 
surface  which  would  compare  well  with  the  similar  definition  based  on  the 
observed  density  structure,  as  is  evident  by  comparing  Figure  25  with  the 
observed  mixed  layer  depth  in  Figure  26. 

2.  Comparison  of  the  Model -Computed  Mixed  Layer  Depth  with  Observed 
Mixed  Layer  Depth 

Figure  26  is  a  time  series  of  the  depth  of  the  model  computed 
turbulent  boundary  layer,  -h(t),  and  of  the  observed  mixed  layer  depth, 
-hg(t) .  The  value  of  h  should  be  less  than  or  equal  to  the  mixed  layer 
depth  as  determined  from  observed  9p/9z  if  the  model  simulation  of  the 

2  Smoothing  was  accomplished  by  giving  the  center  data  point  a  weight 
of  2  and  all  of  the  surrounding  data  points  a  weight  of  1. 
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vertical  density  structure  changes  are  accurate.  Figure  27  is  a  plot  of 
the  same  functions,  h(t)  and  hg(t)  between  time  =  200  hours  and  time  =  400 
hours.  Notice  that  the  model  results  compare  favorably  with  the  observa¬ 
tions  except  for  a  period  of  about  one  day  starting  at  hour  285.  Figure 
28  is  the  model  computed  mixed  layer  depth,  h^(t).  The  value  of  -h-j(t) 
is  taken  to  be  the  depth  at  which  the  model  is  0.02atirits  greater  than 
the  value  of  <ot>  where  <  >  denotes  the  vertical  average  for  the  turbu¬ 
lent  boundary  layer.  This  definition  is  equivalent  to  that  for  hg(t)  if 
the  turbulent  boundary  layer  is  homogeneous.  Figure  29  shows  both  h^(t) 
and  the  observed  mixed  layer  depth,  hQ ( t )  for  the  same  time  period  as 
Figure  28.  Notice  that  h-j£  h  for  all  t.  Although  hg(t)  compares  well 
with  h(t),  there  are  times  during  which  hg(t)  is  more  comparable  to  h-j(t). 
For  example,  for  the  period  between  t  =  285  hours  and  t  *  310  hours, 
hQ(t)  is  most  similar  to  h-| ( t ) .  This  happens  because,  although  the  tur¬ 
bulent  boundary  layer  may  be  quite  shallow,  the  associated  transient 
thermocline  is  not  sufficiently  strong  to  influence  the  values  for  h-j(t) 
and  hQ(t). 

Only  during  one  time  period,  at  about  time  =  100,  is  the  observed 
mixed  layer  depth  found  to  be  significantly  deeper  than  either  the  tur¬ 
bulent  boundary  layer,  h(t),  or  the  mixed  layer  depth,  h-|(t),  computed 
by  the  model.  This  difference  between  the  predicted  and  observed  mixed 
layer  depth  could  be  the  result  of  R/P  FLIP  drifting  through  a  warm  core 
eddy  or  filament  of  water  mass  having  more  tropical  characteristics. 
Figures  21  and  22,  which  are  the  graphs  of  the  advective  components  of 
the  heat  and  salt  budgets,  show  a  large  increase  in  the  magnitudes  of 
both  the  advection  of  heat  and  salt  at  about  time  =  100  hours.  This 
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large  increase  in  advection  of  both  heat  and  salt  is  consistent  with  the 
hypothesis  that  R/P  FLIP  drifted  across  either  an  eddy  or  a  filament  of 
water  from  the  south  of  the  Subtropical  front.  The  Subtropical  water 
mass  which  lies  to  the  south  is  warmer  and  more  saline  than  the  Eastern 
North  Pacific  Central  water. 


3.  Salinity  Effects  on  the  Mixed  Layer  Depth 

An  additional  model  simulation  was  conducted  in  order  to  deter¬ 
mine  what  effect  precipitation  has  on  the  mixed  layer  depth.  A  steady 

-5  -2  -1 

rate  of  precipitation  of  5  x  10  gm  cm  sec  was  prescribed  over  a  12- 
hour  period  starting  at  t  =  200  hours.  At  all  other  times  precipitation 
was  assumed  to  be  negligible.  This  event  corresponds  in  magnitude,  at 
least,  with  the  rate  of  precipitation  reported  by  Simpson  [1977]  for 
5  February  1974.  Figure  30  shows  a  rapid  decrease  in  the  surface  salin¬ 
ity  with  the  injection  of  the  precipitation.  It  drops  from  about  34.11  °/00 
to  33.99  °/00.  This  drop  in  surface  salinity  causes  a  rapid  rise  in  the 
surface  salinity  flux  as  shown  in  Figure  31.  Figure  32  shows  the  model 
prediction  of  the  mixed  layer  depth  with  and  without  the  precipitation 
event.  The  associated  increase  in  downward  buoyancy  flux,  b'w1 ^  due 
to  the  increase  in  salinity  flux  (Figure  31)  caused  the  model  mixed 
layer  depth  to  deepen  much  less  rapidly  than  for  the  case  without  pre¬ 
cipitation.  The  effect  would  have  been  even  more  apparent  if  the  pre¬ 
cipitation  period  had  not  coincided  with  strong  winds. 

Analysis  of  Figure  29  shows  that  the  turbulent  boundary  layer 
depth  predicted  by  the  model  was  approximately  20  meters  deeper  than  the 
observed  mixed  layer  depth  at  t  5  210  hours.  This  time  period  coincides 
with  that  reported  by  Simpson  and  Paulson  [j979^]  as  the  period  of  most 


intense  precipitation.  These  results  suggest  that,  if  the  amount  of 
precipitation  were  accurately  known  and  included  in  the  model  boundary 
conditions,  there  would  be  an  even  better  correlation  between  the  model 
mixed  layer  depth  and  the  observed  mixed  layer  depth. 


III.  SUMMARY 


This  analysis  of  15  days  of  salinity,  temperature  and  density  pro¬ 
files  from  R/P  FLIP  as  part  of  the  NORPAX  POLE  experiment  and  the  ac¬ 
companying  model  study  suggest  that,  non  one-dimensional  effects  were 
very  important  for  both  the  heat  and  salt  budgets.  Vertical  mixing  and 
mixed  layer  depth  are  controlled  primarily  by  the  one-dimensional  response 
of  the  turbulent  kinetic  energy  budget  to  the  local  atmospheric  forcing. 

It  is  apparent  that  the  non  one-dimensional  effects  may  be  largely  ex¬ 
plained  by  the  drift  of  R/P  FLIP.  Results  from  POLE  site  showed  that  this 
apparent  horizontal  temperature  advection  (as  observed  by  the  moving  plat¬ 
form)  was  important  on  all  time  scales  and  dominated  the  small  time  scale. 
Nevertheless,  the  one-dimensional  model  computation  of  mixed  layer  depth 
compared  favorably  with  the  observed  mixed  layer  depth.  It  is  interest¬ 
ing  that  even  though  this  apparent  advection  usually  dominated  the  salt 
budget  over  both  the  long  and  short  time  scales,  the  one-dimensional  mixed 
layer  model  response  to  the  surface  salinity  flux  can  be  significant. 
During  the  single  period  of  strong  precipitation,  the  mixed  layer  depth 
was  shallower  (on  the  order  of  tens  of  meters)  than  it  would  have  been 
without  precipitation.  This  is  consistent  with  the  results  of  Miller 
[1976] ,  showing  that  salinity  flux  can  at  least  occasionally  play  an 
important  part  in  mixed  layer  dynamics  in  the  subtropics.  Further 
studies  would  have  to  be  conducted  with  observational  data  that  included 
the  time,  duration,  and  amount  of  precipitation  before  a  firm  conclusion 
can  be  made  on  the  longer-term  importance  of  salinity  effects  on  the 
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mixed  layer  depth.  Qualitatively,  model  performance  was  improved  when 


the  influence  of  precipitation  was  taken  into  account. 

Simpson  and  Paulson  [l979|  concluded  that  the  one-dimensional  models 
failed  to  predict  the  mixed  layer  depth  and  temperature  observed  during 
POLE.  Why  did  the  models  of  Pollard  et  al  ^973^  and  Niiler  97s|  ,  which 
were  used  by  Simpson  and  Paulson,  fail  to  predict  the  observed  mixed  layer 
depth?_  One  part  of  the  explanation  for  the  model's  apparent  failure  could 
be  the  definition  of  mixed  layer  depth  used  by  Simpson  and  Paulson.  An 
accurate  determination  of  the  mixed  layer  depth  from  a  vertical  profile 
can  only  be  made  by  considering  the  entire  profile,  an  erroneous  analysis 
of  the  mixed  layer  depth  is  probable  whenever  there  is  a  transient 
pycnocline  above  2=5  meters.  Therefore,  contrary  to  the  results  of 
Simpson  and  Paulson  [l979]  ,  it  has  been  shown  that  the  one-dimensional 
turbulence  closure  model  used  here  does  seem  to  explain  most  of  the  ob¬ 
served  mixed  layer  depth  changes  and  that  the  mixed  layer  depth  is  more 
influenced  by  the  local  surface  fluxes  than  by  horizontal  advection. 
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Figure  1.  The  location  of  the  R/P  FLIP  during  the  POLE 

experiment  (35°N,  155°W)  in  relation  to  general 
oceanic  features.  (Simpson  and  Paulson  1977) 
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Figure  2.  Idealized  density  and  mean  velocity  profiles  of  the  ocean 

mixed  layer  ( - ).  Mixed  layer  depth  is  (h),  (6)  is  the 

thickness  of  the  interface  or  entrainment  zone.  Region  I 
is  the  fully  turbulent  mixed  layer  depth.  Region  II  is  the 
slightly  stable,  intermittently-turbulent  entrainment  zone. 
Region  III  is  the  stable  underlying  watermass  having  negli¬ 
gible  vertical  fluxes  in  comparison  to  those  of  Region  I. 
(Garwood  1977) 
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Figure  3.  Typical  POLEX  data  temperature,  salinity,  and  sigma-t  pro¬ 
files  taken  from  Simpson  and  Paulson  [1979].  Added  to  the 
profile  is  Region  I  the  fully-turbulent  mixed  layer  of 
depth  h.  Region  II  the  slightly  stable,  intermittently- 
turbulent  entrainment  zone  of  thickness  6.  Region  III  is 
the  stable  underlying  water  mass  having  negligible  fluxes 
in  comparison  to  those  of  Region  I. 
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Figure  4 
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Typical  P01EX  data  temperature,  salinity,  and  sigma-t  pro¬ 
files  taken  from  Simpson  and  Paulson  [1979J.  Added  to  the 
profile  is  Region  I  the  fully-turbulent  mixed  layer  of 
depth  h.  Region  II  the  slightly  stable,  intermittently- 
turbulent  entrainment  zone  of  thickness  6.  Region  III  is 
the  stable  underlying  water  mass  having  negligible  fluxes 
in  comparison  to  those  of  Region  I. 
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Figure  5.  Mechanical  energy  budget  for  the  ocean  mixed  layer.  Asterisks 
indicate  those  processes  that  must  be  parameterized  to  close 
the  system  of  equations.  (Garwood  1977) 
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computed  hourly  values  of  the  total  downward  heat 
it  the  surface.  The  origin  of  the  abcissa  is  1000 
I  January  1974. 
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Figure  17.  Interpolated  hourly  values  of  mass  and  buoyancy  fluxes, 
surface  salinity,  sea  surface  temperature,  surface 
sigma-t,  and  the  various  components  of  the  heat  flux 
involved  in  the  determination  of  the  buoyancy  flux. 
(Simpson  and  Paulson  1977) 
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Figure  18.  Interpolated  hourly  values  of  the  various  components 
of  the  surface  momentum  and  heat  balance  are  shown. 
(Simpson  and  Paulson  1977) 
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Figure  19.  The  observed  changes  in  heat  content  in  the  upper 
40  meters  (solid  line).  Model  computed  accumulated 
surface  heat  flux  (dotted  line).  The  origin  of  the 
abscissa  is  0000  hour  (local)  30  January  74. 
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Figure  20.  The  observed  changes  in  salt  content  in  upper  40 
meters  (solid  line).  Model -computed  accumulated 
surface  salinity  flux  (3TwT),  dashed  line,  which 
is  negligible  in  comparison  with  the  observed 
changes  in  salt  content.  The  origin  of  the 
abscissa  is  0000  hour  (local)  30  January  74. 
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Figure  21.  Advection  of  heat  (solid  line)  determined  by 
subtracting  the  accumulated  heat  flux  from  the 
change  in  heat  content  (Figure  21).  Least 
squares  fit  to  data  resulted  in  a  slope  2.0  cal 
cm'^  hour’  . 
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Figure  22.  Advection  of  salt  (solid  lire)  determined  by 
subtracting  the  accumulated  surface  salinity 
flux  from  the  changes  in  salt  content  (Figure 
20).  Least  square  fit  to  data  resulted  in  a 
slope  0.45  0/00  cm  hour-^. 
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Figure  23.  Comparison  of  model  computed  mixed  layer  temperature  (solid  line)  with 

observed  for  sea  surface  temperature  (dotted  line).  Model  computed  mixed 
layer  temperature  has  been  offset  slightly  for  the  comparison.  The  origin 
of  the  abscissa  is  1200  GMT  30  January  1974. 
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Figure  24.  Isotherm  distribution  for  the  period  30  Jan  74  through 
14  Feb  74.  The  contour  interval  is  0.5°C. 
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Figure  26.  The  dashed  line  with  the  symbol  (0)  is  the  observed  mixed  layer  depth 
hg(t).  The  solid  line  is  the  model  computed  turbulent  boundary  layer 
depth,  h(t).  The  origin  of  the  abscissa  is  1000  GMT  28  January  1974. 


The  dashed  line  with  the  symbol  (0)  is  the  observed  mixed  layer 
depth,  hQ(t).  The  solid  line  is  model  computed  turbulent 
boundary  tayer  depth,  h(t). 
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Figure  29.  Model  computed  mixed  layer  depth  (solid  line),  h.|(t).  Observed 
mixed  layer  depth  (dashed  line),  hg(t). 
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Figure  31.  Model  computed  hourly  values  of 
with  hypothetical  precipitation 
starting  at  time  =  220. 
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